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Abstract 

The dimensionally regularized master planar double box Feynman dia- 
gram with four massive and three massless lines, powers of propagators 
equal to one, all four legs on the mass shell, i.e. with p 2 = m 2 , i = 
1,2,3,4, is analytically evaluated for general values of m 2 and the Man- 
delstam variables s and t. An explicit result is expressed in terms of 
polylogarithms, up to the third order, depending on special combina- 
tions of m 2 , s and t. 
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1. Systematical analytical evaluation of two-loop Feynman diagrams with four 
external lines within dimensional regularization |3J began two and a half years ago. 
In the pure massless case with all end-points on-shell, i.e. p\ = 0, i — 1,2,3,4, 
the problem of analytical evaluation of two-loop four-point diagrams in expansion 
in e = (4 — d)/2, where d is the space-time dimension, has been completely solved 
in 0, |3|, U ||, ||]. Any such diagram can be expressed, in Laurent expansion in e 
up to a finite part, through polylogarithms and generalized polylogarithms up to 
the fourth order, depending on the ratio of the Mandelstam variables s and t. The 
corresponding analytical algorithms have been successfully applied to the evaluation 
of two-loop virtual corrections to various scattering processes in the zero-mass 
approximation. 

If in addition to the Mandelstam variables s and t there is one more massive 
parameter, the four-point diagrams become much more complicated. However, in 
the case of one leg off-shell relevant to the process e + e~ — > 3jets (see, e.g., H), the 
problem of the evaluation has been solved: master integrals have been evaluated 
H pD| , either by use of Mellin-Barnes representation or the method of differential 
equations fLl"] , and a reduction procedure has been developed fLO] . (See |12| where 
the present status of NNLO calculations of the process e + e~ — > 3jets is characterized.) 

The purpose of this paper is to turn attention to on-shell four-point diagrams 
with a non-zero internal mass and analytically evaluate the scalar double box diagram 
shown in Fig. 1. The calculational experience, in particular obtained in the above 
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Figure 1: Planar double box diagram. Solid and dashed lines denote massive and 
massless propagators, respectively. 

mentioned works on the evaluation of four-point diagrams, tell us that if such master 
integrals can be evaluated, the problem can be also completely solved, after evaluating 
other master integrals and constructing a recursive procedure that expresses any given 
Feynman integral with general numerators and integer powers of propagators through 
the master integrals. Therefore this explicit analytical result can be considered as a 
kind of existence theorem, in the sense that it strongly indicates the possibility to 
analytically compute various scattering processes in two loops without putting masses 
to zero. 
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To arrive at this result we derive in the next section an appropriate Mellin-Barnes 
(MB) representation for the general planar double box within dimensional regulariza- 
tion. Then we turn to the master double box and use a standard procedure of taking 
residues and shifting contours to resolve the structure of singularities in the parameter 
of dimensional regularization, e. The initial MB integral is eventually decomposed 
into seven pieces where expansion of the integrand in e becomes possible. After eval- 
uating these expanded integrals we obtain an explicit analytical result expressed in 
terms of polylogarithms, up to the third order, depending on special combinations of 
m 2 , s and t. In the last section, leading order terms of the asymptotic expansion of the 
considered diagram in the limit of small m obtained by the strategy of expansion by 



regions [13, 14, 15] are presented. They serve as a crucial check of the result obtained. 



We conclude with a discussion of the results of the paper and open problems. 

2. Let us consider the general on-shell double box diagram of Fig. 1, i.e. with 
general irreducible numerator and powers of propagators. We choose this irreducible 
numerator and the routing of the external momenta as in |J. For convenience, we 
consider the factor with (k + pi + p 2 +P3) 2 corresponding to the irreducible numerator 
as an extra propagator but, really, we are interested only in the non-positive integer 
values of a 8 . This general double box Feynman integral takes the form 

d d k d d l 



B(ax, . . . , as', s, t, m 



2. 



(k 2 - m 2 ) a i [(k + Pi) 2 )} a2 [{k + pi+ P2) 2 - m 2 }^ 



[jk + Pi + p 2 + p 3 ) 2 ]" as m 

X [(l+p 1 +p 2 ) 2 -m 2 )}^[(l + p 1 +p 2 +p 3 ) 2 }^(l 2 -m 2 )^[(k-l) 2 }^ ' 1 ' 

where s = {p\ + P2) 2 , t = (P2+P3) 2 , and k and I are respectively loop momenta of the 
left and the right box. Usual prescriptions k 2 = k 2 + iO, s = s + iO, etc. are implied. 

To resolve the singularity structure of Feynman integrals in e it is very useful to 
apply the MB representation 

1 1 r +ioo yz 

dz— -Y{v + z)Y{-z) , (2) 



(X + Y) v T{v) 2m J-ioo X v + Z 

that makes it possible to replace sums of terms raised to some power by their products 
in some powers, at the cost of introducing an extra integration. In 0, |4], |9| MB 
integrations were introduced directly in alpha/Feynman parametric integrals. It turns 
out more convenient to follow (as in 0) the strategy of [[TEJ and introduce, in a suitable 
way, MB integrations, first, after integration over one of the loop momenta, I, and 
complete this procedure after integration over the second loop momentum, k. In 
fact, the procedure of is straightforwardly generalized by introducing two extra 
MB integrations when separating terms with m 2 after each of the integrations over 
the loop momenta, and after appropriate changes of variables leads to the following 
sixfold MB representation of ([!]): 

^/ 2 ) 2 (-ir 



B(ai, . . . , a 8 ; s,t,m 
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{2-Kif J-ioo W Z3 \s) \s) T( ai + z 3 + z 4 )T(a 3 + z 2 + z 4 ) 
T(4 - oi3 - 2a 28 - 2e + z 2 + Z3)r(ai 238 - 2 + e + z 4 + z 5 )T(a 7 + w - z 4 ) 

X ~ " 

T(4 — a 4 6 — 2a 57 — 2e — 2w — 2z\ — z 2 — 23) 
r(a 4 5 67 - 2 + e + u; + 2;i- z i )T(a s — z 2 — z 3 — z A )Y(-w - z 2 - z 3 - z 4 ) 

r(4 - ai238 - 2e + 10 - z 4 )r(a 8 - 10 - z 2 - z 3 - z 4 ) 
T(a 5 + w + z 2 + z 3 + %)r(2 - a 567 - e - w - 21 - z 2 ) 
r(4 - ai 3 - 2a 28 - 2e + 2 2 + 23 - 2z 5 ) 
xT(2 - a 457 - e - w - z x - z- 3 )Y{2 - a X2 % - e + z 2 - z^)Y{2 - a 238 - e + z 3 - z 5 ) 
xr(4 - a 46 - 2a 57 - 2e - 2w - z 2 - z 3 )T{-z x )T{-z h ) , (3) 

where a 4 56 7 = a 4 + a 5 + + a 7 , a i3 = ai + a3, etc., and integration contours are chosen 
in the standard way. 

In the case of the master double box, we set a$ = 1 for i — 1, 2, . . . , 7 and a s = 
and obtain 

B ( °)(s,t,m 2 ;e) = 5(1, . . . , 1, 0; s, t, m 2 ; e) 

H /2 ) 2 1 r ri A dz /m^V 1+z 7ty r(l + «,)r(-«,) 

r(-2e)(-s) 3 + 2e (27ri) 6 i-ioo / = \ "V" 5 / ^ T(l - 2e + w; - ^ 4 ) 
T(2 + e + w + Z\ — z 4 )r(— 1 — e — w — Z\ — z 2 )T(—l — e — w — Z\ — z 3 )T(—Zi) 

r(i + z 2 + z 4 )r(i + -2 3 + ^ 4 )r(-2e + z 2 + ^ 3 - 225) 

^ T(-e + z 2 - z 5 )T(-e + z 3 - z 5 )r(l + e + z A + z 5 )r(-^ 5 )r(-2e + ^ + Z3) 

X T(-2 -2e-2w- 2z x - z 2 - z 3 ) 

xT(-2 - 2e - 2w - z 2 - z 3 )Y{\ + w + z 2 + z 3 + z i )T{z 2 + z 4 )T(z 3 + z 4 ) 

xT(l + w - z i )T(-z 2 - z 3 -z A ). (4) 

Observe that, because of the presence of the factor r(— 2e) in the denominator, we 
are forced to take some residue in order to arrive at a non-zero result at e = 0, so 
that the integral is effectively fivefold. 

The resolution of singularities in e is performed also in the standard way (see 
0. f|, D, §]) and reduces to shifting contours and taking residues. The goal of this 
procedure is to decompose a given integral into pieces where the Laurent expansion e 
of the integrand becomes possible. This is how such procedure can be performed for 

@: 

1. Take minus residue at z 3 = —2 — 2e — 2w — z 2 , then minus residue at w = — 1 — 2e, 
then residue at z 4 = 0, then residue at z 2 = 0, expand in a Laurent series in e 
up to a finite part. Let us denote the resulting integral over Z\ and z 5 by B x . 

2. Take minus residue at z 3 = —2 — 2e — 2w — z 2 , then minus residue at w — — 1 — 2e, 
then residue at z 4 = 0, and change the nature of the first pole of T(z 2 ) (choose 
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a contour from the opposite side, i.e. the pole z 2 will be now to the right of the 
contour), then expand in e. Denote this integral over zi, z 2 and z§ by B 2 . 

3. Take minus residue at z 3 = —2 — 2e — 2w — z 2 , then minus residue at w = — 1— 2e, 
then change the nature of the first pole of T(z4), then take a residue at z 2 = —£4, 
then take a residue at Z4 = — e and expand in e. This resulting integral over z\ 
and z 5 is denoted by B 3 . 

4. Take minus residue at z 3 = —2 — 2e — 2w — z 2 , then minus residue at w = — 1 — 2e, 
then change the nature of the first pole of r(^ 4 ), then take a residue at z 2 = — z±, 
then change the nature of the first pole of T(2(e + 214)) and expand in e. The 
resulting integral over z\, z 4 and z§ is denoted by B4. 

5. Take minus residue at z 3 = —2 — 2e — 2w — z 2 , then minus residue at w = — 1 — 2e, 
then change the nature of the first pole of T(z 4 ), then change the nature of the 
first pole of T(z 2 + 24) and expand in e. The resulting integral over z 1; z 2 , z 4 and 
Z5 is denoted by B 5 . 

6. Take minus residue at z 3 = —2 — 2e — 2w — z 2 , then change the nature of the 
first pole of T(— 2(1 + 2e + w)), then take minus residue at z 4 = 1 + w, then 
minus residue at z 2 = — 1 — 2e — w and expand in e. The resulting integral over 
w, zi and z 5 is denoted by B 6 . 

7. Change the nature of the first pole of T(— 2 — 2e — 2w — z 2 — z 3 ), then take minus 
residue at Z4 = —z 2 — z 3 , then a residue at z% = 2e — z 2 , then take a residue at 
z 2 = 2e and expand in e. The resulting integral over w, z\ and z^ is denoted by 
B 7 . 

One can see that all other contributions vanish at e = 0. By a suitable change 
of variables, one can observe that Bj = Bq. In fact, the dependence of the first five 
contributions on the Mandelstam variable t is trivial: they are just proportional to 

1/*- 

The two-dimensional integrals Bi and B 3 are products of one- dimensional integrals 
which are taken by closing contour to the left and summing up resulting series with 



the help of formulae related to those of Appendix B.2 of ||17|| . The three-dimensional 
integral B 2 is evaluated by closing the integration contours over z\ and z 5 to the left, 
summing up resulting series and applying a similar procedure to a final integral in z 2 . 
The corresponding result is naturally expressed through polylogarithms, up to Li 3 (a;) 
depending on s and m 2 in terms of the variable 



yjl - s/(4m 2 ) + y/=s/(2m) 
y/l - s/(4m 2 ) - y/=a/(2m) 

The form of this result provides a hint about possible functional dependence of 
results for the three-(four-) dimensional integrals B4 (-B5), and a heuristic procedure 
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which was explicitly formulated in |18| turns out to be successfully applicable here. 
First, all the contributions, in particular i? 4 and B 5 , are analytic functions of s in a 
vicinity of the origin. One can observe that any given term of the Taylor expansion 
can be straightforwardly evaluated because the corresponding integrals over z 2 and 
Z4 are recursively taken. It is, therefore, possible to evaluate enough first terms (say, 
30) of this Taylor expansion. Then one takes into account the type of the functional 
dependence mentioned above, turns to a new Taylor series in terms of the variable 
v — 1 and assumes that the n-th term of this Taylor series is a linear combination, 
with unknown coefficients, of the following quantities: 



n 



1_ S^n)_ S?(n) S 3 (n) 
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where Sk(n) = YHj=i3~ ■ (Here some terms of the fourth order are listed. See [IS 
for more details.) Using information about the first terms of the Taylor series one 
solves a system of linear equations, finds those unknown coefficients and checks this 
solution with the help of the next Taylor coefficients. 

This experimental mathematics has turned out to be quite successful for the eval- 
uation of B 4 and B 5 . Finally, the contribution B e is a product of a one-dimensional 
integral in zi, which is easily evaluated, and a two-dimensional integral in w and Z5 
which involves a non-trivial dependence on t and is evaluated by closing the inte- 
gration contours in z§ to the left, summing up a resulting series in terms of Gauss 
hypergeometric function for which one can apply a parametric representation. After 
that the internal integral over w is taken by the same procedure and, finally, one takes 
the parametric integral. 

The final result takes the following form: 



B^(s,t,m 2 -e) 
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l+2e 
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b 2 (x) 
h{x) 



b i{x) 



where x = 1/Jl — Am 2 /s, y = 1/Jl — Am 2 /t, and 
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-(tt 2 /3)(4/ 2 -m x - 3p x ) + (8/3)Z| + 14C 3 , 

, . /I — x 



(m x - p x ) 
2x 
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-(4/3) [tt 2 - 6Z 2 + 3m 2 + &m x (2l 2 - 2l x - p x ) + 12/^ - 3p 2 ](Li 2 (x) - Li 2 (-x)) 



■(8/3)[7r 2 - 6/ 2 + 6/^ - 6m x (l x + p x - 2/ 2 )]Li 2 
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+8(m x - p x 
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2x 
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+ 2(l x -m x + / 2 )Li 2 



-2x 



1 — x 



-8{m x -p x )(2l x -p x - 5m x + 4l 2 )(-m x p x + l 2 {m x + p x ) - l\ + tt 2 /6) 
-(20/3)m x + (164/3)m^ - 40m 2 p 2 - (A/3)m x p 3 x _ (8/3 ) p 4 
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-Al 2 (7ml + 21m 2 x p x - Am x l xPx - 23m x p 2 x + M x p 2 x - p 3 x ) 
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The last piece of the finite part comes from B 6 and B 7 : 
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(9) 



Here Li a (z) is a polylogarithm [T^]. The following abbreviations are also used: ( 3 = 
C(3), l z — hxz for z = x,y, 2, p z = ln(l + z) and m z = ln(l — z) for z — x,y, xy. 

This result is presented in such a way that it is manifestly real at small negative 
values of s and t. From this Euclidean domain, it can be easily analytically continued 
to any other domain. 

3. The result is in agreement with the leading power behaviour in the 

(Sudakov) limit of the fixed-angle scattering, m 2 -C \s\, \t\. This asymptotics is ob- 
tained by use of the strategy of expansion by regions |L3l |15|. The structure of 



regions is very rich. The following family of seventeen regions participates here: 

(h-h), (lc-h), . . . , (4c-h), (lc-lc), . . . , (4c-4c), 
(lc-3c), (2c-4c), (lc-4c), (2c-3c), 
(luc-2c), (2uc-lc), (3uc-4c), (4uc-3c). 

Here h denotes hard, c - collinear and uc - ultracollinear regions for the two loop mo- 
menta. (See [14| and Chapter 8 of [|15| for definitions of these regions.) In particular, 
the (h-h) contribution is nothing but the massless on-shell double box 0. The con- 
tributions (lc-lc), (3c-3c), (lc-3c) as well as the symmetrical contributions (2c-2c), 
(4c-4c) and (2c-4c) are not individually regularized by dimensional regularization. 
The poles in the auxiliary analytic regularization turn out to be of the second order 
and are cancelled in the sum. After adding the rest of the contributions, the poles of 
the third and fourth order in e are cancelled. Following this procedure, we obtain 



S (0) (s,t,m 2 ;e) 



/ i7r d/2 e -7Eey 




'(2/3)L 3 + (n 2 /3)L + 2( 3 



1 



+ 



-(2/3)L 4 + 2 \n(t/s)L 3 - 2(\n 2 (t/s) + 4tt 2 /3)L 2 
4Li 3 (-t/s) - 4 ln(t/s)Li 2 {-t/s) + (2/3) ln 3 (t/s) - 2 ln(l + t/s) \n 2 (t/s) 



+ (8tt 2 /3) \n(t/s) - 2tt 2 ln(l + t/s) + lOCs] L + vr 4 /36} + 0(m 2 L 3 , e) , (10) 

where L = ln(— m 2 /s). This asymptotic behaviour is reproduced when one starts 
from result 
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The analytical result presented above agrees also with results based on numerical 
integration in the space of alpha parameters |2(J (where the 1% accuracy for the 1/e 
and e° parts is guaranteed). 

Let us stress that, in the present case with a non-zero mass, there are no collinear 
divergences and the poles in e are only up to the second order, so that the resolution 
of singularities in e in the MB integrals is relatively simple. Therefore, it looks 
promising to use the technique presented, starting from @, for the evaluation of any 
given master integral. The construction of a recursive algorithm that would express 
any given planar double box through some family of master integrals is one of the 
next problems, as well as similar problems for the evaluation of massive on-shell non- 
planar double-box Feynman integrals. Another possible scenario in the situation, 
where the mass is small but still not negligible, is to evaluate the leading power (and 
all logarithms) asymptotics when m —>■ 0. Anyway, the (h-h) contribution to any 
such asymptotic behaviour is obtained by the algorithms of 0, |3|, [|, [5|. |(| . 

It is interesting to note that, in the above result, there are no so-called two- 
dimensional harmonic polylogarithms EI[ which have turned out to be adequate 



functions to express results for the double boxes with one leg off-shell fllP] . It is also 
an open question whether this phenomenon takes place for general massive on-shell 
double boxes. 
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